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ABSTRACT 

We present here an alternative method for calculating magnifications in gravitational 
lensing calculations - the Ray Bundle method. We provide a detailed comparison be- 
tween the distribution of magnifications obtained compared with analytic results and 
conventional ray-shooting methods. The Ray Bundle method provides high accuracy 
in the weak lensing limit, and is computationally much faster than (non- hierarchical) 
ray shooting methods to a comparable accuracy. 

The Ray Bundle method is a powerful and efficient technique with which to study 
gravitational lensing within realistic cosmological models, particularly in the weak 
lensing limit. 

Key words: gravitational lensing - methods: numerical 



1 INTRODUCTION 

Gravitational lensing is the study of the effects of matter 
on the propagation of light. The most obvious observational 
results are the production of multiple images (as first seen 
with the multiply i maged quasar 0957+561 by Walsh, Car- 
swell & Weymann ( 197E ) ) , the creation of giant luminous 
arcs (first identified in the ga laxy clusters Abell 370 and 
CI 2244 by Lynds & Petrosi an (1986 ) and, independently, in 
Abell 370 by Soucail et al. (1987)) and the large magnifica- 
tions of source flux seen in microlensing events (for example, 
brightening of a single image of the multiply-imaged quasar 
2237+0305 due to compact objec ts in the lensing galaxy, 
first detected by Irwin et al. (198£)). 



The paths of light rays from a source to the observer are 
conveniently described with the gravitational lens equation 
(Equation (|l]) below). This equation is highly non-linear, so 
that, except for a small number of specific cases, there are no 
analytic solutions. In particular, there is no straightforward 
result which determines the image locations or magnifica- 
tions for an ensemble of many lenses. This presents a serious 
problem when we wish to study the lensing properties of a 
complex lensing structure such as the Universe. 

Fortunately, a number of approximate numerical meth- 
ods have been developed which allow us to calculate magni- 
fications and other properties of a collection of lenses. Fore- 
most amongst these are th e Ray Shooting methods, intro- 
duced by Paczyriski ( 1986 ) and Kayser, Ref sdal fc Stabell 
(1986) and d eveloped by Schnei der fc Weiss (1987), Kayser 
et al. (1989) and Lewis et al. (1993). For a discussion of 
other methods, such as the use of a scalar deflection poten- 



tial, and the o ptical scalar equations, see Schneider, Ehlers 
& Falco ( |l992j ). 

In this paper, we present an alternative technique for 
calculating the magnification properties of an ensemble of 
lenses - the Ray Bundle method (RBM) . The RBM is partic- 
ularly well suited to studies of the weak lensing limit, where 
we are not concerned with the creation of multiple (com- 
parably bright) images. Like the Ray Shooting method, the 
RBM uses backwards propagation of light rays from the ob- 
server to the source, which are deflected by the distribution 
of lenses, and are mapped to the source plane. Whereas the 
Ray Shooting method collects the deflected light rays within 
a rectangular grid of pixels, we consider an infinitesimal bun- 
dle of rays which form a circular image, and maintain this 
association to the source plane. 

At first sight, the Ray Bundle method may seem to 
be less computationally efficient than the Ray Shooting 
method, as a grid bas ed technique lends i tself to (fast) hier- 
archical calculations (Wambsganss 1990). However instead 
of requiring > 100 light rays per source grid to reduce the 
statistical error, we need only use a bundle of N lay = 8 rays 
to obtain magnifications which are correct to better than 5 
per cent for an equivalent source size. 

With a Ray Shooting method, we have both image and 
source pixels, but cannot easily determine the correspon- 
dence between them. By keeping track of the individual light 
bundles with the RBM, we are able to monitor the shape dis- 
tortions of the beam caused by the shear and convergence 
of a lens ensemble. This is of particular interest when the 
RBM is applied to multiple lens plane geometries (as are 
conventionally used for studying cosmologically distributed 
lenses). Details of the beam shape allows for the opportu- 
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nity to make comparisons between results calculated with 
the gravitational lens equation, and those using the optical 
scalar equations (which are more easily applied to smooth 
mass distributions, or appr oximate mass distributi ons such 
as the Swiss cheese model (Einstein & Straus 1945)). 

In Section ^ we introduce various basic results of grav- 
itational lensing, particularly with regards to magnification 
and the magnification probability distribution. The Ray 
Bundle method is introduced in Section^ and compared 
in detail with analytic solutions for the Schwarzschild lens 
model. By obtaining the magnification probability distribu- 
tion with both the Ray Bundle and Ray Shooting methods 
for a variety of lens geometries, we demonstrate the general 
applicability of the RBM. We do not discuss applications 
of the RBM here, but reserve details for Fluke, Webster & 
Mortlock (in preparation) where the Ray Bundle method is 
used to investigate weak lensing within realistic cosmological 
models (generated with N-body simulations). 
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Figure 1. Geometrical arrangement for the gravitational lens 
equation for an observer at O. £ is the impact parameter in the 
lens plane, r) the sources position, and &(£) the deflection angle of 
the light ray. The distances (-Dij) are angular diameter distances 
between the [ojbserver, [djeflector and [sjource planes. 



2 GRAVITATIONAL LENSING 

We present here a number of important results from gravi- 
tational lensing which we will require: the gravitational lens 
equation, the magnification of source flux and the magnifica- 
tion probability distribution. Several excellent sources exist 
which describe gravitational lensing far more comprehen- 
sively than may be discussed here, for e xampl e Schneider 
et al. ( |l992| ) and Narayan & Bartelmann (^99§). 



2.1 The gravitational lens equation 

The deflection of a light ray by a massive object is con- 
veniently expressed with the gravitational lens equation 
(GLE), which may be derived from simple geometrical ar- 
guments (as demonstrated in Figure]]]). 

The GLE relates the impact parameter, £, in the lens 
(deflector) plane, the source position, rj, in the source plane 
and the deflection angle, a(£), of the light ray: 

L>, ■ 



od 



(1) 



The distances (£^ij) in Equation (gj) are angular diameter 
distances between the [o]bserver, [djeflector and [sjource 
planes. For a given source position, an image will occur for 
each value of £ which is a solution of Equation (|l|) . We de- 
fine the lens axis as the line from the observer through a 
single lens (at the origin of the lens plane), and which is 
perpendicular to both the lens and source planes. 

The GLE may be recast into a dimensionless form by 
introducing the scaling lengths £o and 770 = ^oDos/D^. 
Defining x — £/fo, D = v/Vo and the dimensionless deflec- 
tion angle 



a(x) 



-Dos Co 
it follows that 
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2.2 Magnification 

For a bundle of light rays passing through a transparent lens, 
the number of photons is conserved, ie. gravitational lensing 
does not change the specific intensity of the source. A change 
in flux, however, can occur as the cross-sectional area of a 
bundle of light rays will be affected by a gravitational lens. 
The change in the apparent luminosity is entirely due to 
the change in the solid angle that the image covers (at the 
expense of the rest of the sky) - with a lens present, more 
photon trajectories are brought to the ob server's eye than if 
there were no lens (Dyer & Roeder 1981). 

If the flux at a frequency v is S v = 7„df2 bs, where I v is 
the specific intensity and df2 b s is the solid angle subtended 
by the source at the observer's location, then the magnifi- 
cation is 



S„ 



IMI SI dfii 



(4) 



with primes denoting quantities when the lens is absent. 

It is usual to measure the magnification with respect 
to either a 'full beam' or an 'empty beam'. The full beam 
refers to the case where matter is smoothly distributed ev- 
erywhere (inside and outside of a bundle of light rays), so 
that there is no shear and the magnification is due entirely to 
convergence. An empty beam corresponds to the case where 
there is a smooth distribution of matter (on average) exter- 
nal to the beam, and no matter within the beam. The empty 
beam is the maximally divergent beam, as there is now no 
magnification due to convergence, and the minimum mag- 
nification is /iempty = 1. The addition of material outside 
of the be am will result i n a total magnification /i > 1 due 
to shear (Schneider 1984). In the work that follows, we will 



calculate magnifications with respect to an empty beam. 



2.3 The magnification probability distribution 

The (differential) magnification probability, p(fi, z)dfi, is the 
probability that the total magnification of a source at red- 
shift z will lie in the range /i to + d/i. The probability is 
subject to the constraints 
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p{n, z)djj, = 1 



j i 
and 



p(n,z)fj,dn = (n)(z) 



(•>) 



(6) 



where (p)(z) is the mean magnification. For a universe where 
the matter is smoothly distributed, (p)(z) = 1, while for 
the empty beam magnifications considered here, {/-i)(z) > 1. 
Equation (||) provides the normalisation of the prob ability, 
while Equation Q expresses the conservation of flux ( Wein- 
berg 19||). 

These constraints are only strictly true when we con- 
sider the magnification probability over the whole sky - in 
the tests we conduct here, we use artificial lens distributions 
that do not cover the entire celestial sphere. 

The probability distribution provides a straightforward 
way of making a comparison between different techniques for 
solving the GLE. Although we might expect the calculation 
of the magnification along a particular line-of-sight to vary 
slightly (depending on source geometry, etc), the statistical 
properties of a distribution of sources should be essentially 
independent of the method. 

In practice, we will use a histogram of the probability as 
a function of fi. This involves summing over a range of mag- 
nifications from to jj.2 = fii + A/i, so that each histogram 
bin represents 



p(fXi,z)Afi = / p(p,z)dfi 



(T) 



In the following, we will refer to this as the magnification 
probability histogram (MPH). 



2.4 Solving the gravitational lens equation 

The GLE is highly non-linear, and in general there are mul- 
tiple solutions for the image locations for a single source 
position. As a result, analytic solutions (where we can in- 
vert the GLE to solve for all £ as a function of 77) exist 
for only special lens geometries and models, which may not 
always be realistic or practical. Instead, various numerical 
methods have been developed which enable solutions to the 
lens equation to be found accurately and efficiently. 

Some of the most widely used methods are those based 
on the concept of Ray Shooting. The basic principle of the 
Ray Shooting method (RSM) is to use Equation (jjj) to prop- 
agate light rays backwards from the observer through a se- 
quence of one or more lens planes to the source plane, where 
the rays are collected on a rectangular pixel grid. Typically 
~ 10 6 pixels are required, with an average of JV ~ 100 rays 
per source plane pixel. The magnification in each source 
pixel is then proportional to the density of rays col- 

lected therein: 



Mpixel(i,j) = A / 'collcctcd(i, j)/N. 



(8) 



A shooting region in the lens plane is chosen such that 
only a few rays from outside of this grid are mapped onto 
one of the sources, otherwise there would be missing flux. 
By using a uniform gri d of image rays, th e relative error is 
approximately N~ 3 ^ 4 (Kayser et al. 1986), which is better 
than the Poisson error if a random distribution of image 



rays was used (however a regular grid may introduce sys- 
tematic errors). As a consequence of this error, it is possible 
to get magnifications /i p i xc i(*, j) < 1 , vio lating the condition 
on the total magnification (Section 2.2), which is entirely a 
numerical effect. 

Early versio ns of this method were introduced by 



Paczynski (1986) and Kayser et al. (1986) and developed 



by Schneider & Weiss (1987) and Kayser et al. (1989). Hier- 



archical tree met hods w ere applied to microlensing scenarios 
by Wa mbsganss ( |l990| ) and Wambsganss, Paczynski & Katz 
( [l990l ). The hierarchical methods approximate the effects of 
lenses which are far from a light ray, and allow the inclusion 
of many thousands of lenses at a low computational cost 
(0(7Vlog 2 iV) for a tree-code with N lenses, versus 0(N 2 ) 
when the contribution of every lens is explicitly calculated). 
An improvement on conventional Ray Shooting methods for 
obtaining statistical properties of microlensing light curves 
can be made with the efficient one-dimensional contour fol- 



lowing algorithm of Lewis et al. (1993). We now introduce 
the Ray Bundle method as an alternative to the RSM. 



3 THE RAY BUNDLE METHOD 

The Ray Bundle method (RBM) is similar to the RSM, in 
that the lens equation is used to propagate light rays back- 
wards from the observer to the source plane. However we 
now consider a bundle consisting of a central ray (the null 
geodesic) surrounded by N ray light rays, which create an 
image shape (usually circular). As the ray bundle passes 
through the lens distribution, its shape will be distorted due 
to shear (stretching along an axis) and convergence (focusing 
due to matter within the beam). For an 'infinitesimal' ray 
bundle, the magnification is determined from Equation (U|) 
by calculating the area of the bundle in the image plane 
(d£l' obs ) and the source plane (dQ obs ). 

Since we are using backwards ray-tracing through a sin- 
gle image position, we do not know where other images may 
occur - hence any measurement of the magnification using 
the RBM will underestimate the total magnification. How- 
ever, this is only a significant problem when we consider im- 
ages located near the critical curves, when the contribution 
to the total magnification due to any other images becomes 
important (see below). The RBM was developed for appli- 
cations in the weak lensing limit, and should be used with 
caution for strong lensing cases. 

3.1 Comparison with analytic solutions 

We can investigate the validity of the RBM by compari- 
son with the various analytic solutions which exist for the 
Schwarzschild lens (see Appendix ^|for a summary). 

Consider first a circular source of radius i? a with centre 
at y c — (yi, c, 2/2, c). The circumference of the source is then 
described by the set of vectors y = (2/1,2/2) with 



yi = 

2/2 = 



2/i,c + R s cos(fli) 
2/2,c + Rs sin(^) 



(9) 



where < cj> < 2tt. 

For each y , w e can solve for the two solutions, x±, 
with Equation (A3). In this case we are using the GLE to 
map from the source plane to the image plane. A source far 



© 1998 RAS, MNRAS 000, §-| 



4 C.J. Fluke et al. 




Analytic source and images 
Ray bundle image and source 
Einstein ring 



-1 1 

Figure 2. Shapes and locations of a circular source, for which 
the boundary points may be determined for both images (solid 
line), and a circular image and corresponding source (short dashed 
line). The lens is located at the cross, and the long dashed line 
is the Einstein radius. In both cases (solid and short dashed line) 
the source is the shape which intersects the Einstein radius. The 
scale is in units of the Einstein radius. 



from the lens axis produces one highly demagnified image 
(/ifaint) located near the lens axis (at a;f a i nt ). The second 
image will have a magnification /^bright > 1, and an angu- 
lar position near the source (at alright)- As the source is 
moved towards the lens axis, the images are stretched in the 
tangential direction^ and become comparable in brightness 
(iMfatatl ~ iMbrightl > 1 as y c — > 0). When y c = 0, the 
two images merge into a highly magnified ring (th e Ein stein 
ring) with total magnification given by Equation (A8). 

Now consider a circular image of radius Ri centred on 
the location of the bright image, x c — alright (y c ), with cir- 
cumferential points x — (xi,X2) 



Xl 
X2 



Xl,c + Ri COs((/>) 

X2, c + -Ri sin(0). 



(10) 



This time, the GLE maps in the opposite direction - from 
the image plane to the source plane. The source shape we 
obtain is stretched along the radial direction, and differen- 
tially compressed in the tangential direction. 

Figure ^| demonstrates the differences between the 
shape and locations of a circular source (solid line), for 
which the circumferential points may be determined for both 
images, and a circular image and its corresponding source 
(short dashed line). In this example, we have considered a 
source which is near the Einstein radius where strong lensing 
effects dominate, and there may be a significant contribu- 
tion to the flux from the second image. In all cases where 
we will apply the RBM, the image is chosen to be well away 
from the Einstein radius (critical curve) , and so the flux lost 
from the second image is not important, as we now show. 



* The Schwarzschild lens has a single (degenerate) caustic point 
at y = 0, which corresponds to a tangential critical curve at x = 1 
(the Einstein ring). 



3.2 The magnification deficit 

We now look at how accurately the RBM approximates the 
total magnification, even though it includes the contribution 
of only one image. If we set R s — Ri, as we expect only small 
changes to the shape and hence radius of the image in the 
weak lensing limit (|a5 c | ^ 1), then this is a two parameter 
problem (Ri,N lay ). 

Defining the RBM magnification in terms of the ray 
bundle image and source areas (A^rbm, j4 s ,rbm) as 



^RBM = 



and the true (total) magnification as 

-<4faint 4* -^bright 



p-true — 



(11) 



(12) 



where j4f a i nt , ^bright are the areas of the two images, then 
the relative error in /jrbm is 

A/iRBM _ \fJ-txae — MB.BM | 

/Atrue /^truc 

Due to the circular symmetry of the Schwarzschild lens 
model, we need only determine the radius, x cut , within which 
the RBM produces a relative error Amrbm > p per cent. 

r MRBM 1 v 

By using iVr ay rays in the image and source bundles, 
we are approximating the shape of a circular image/source 
by a polygon with N lay sides. Clearly, when N lay 3> 1, we 
will have a reasonable approximation to the true shape of 
the image/source. However, to improve the speed of the ray 
bundle method (at the cost of a small error), we ideally 
want to select a small value of iV ray ( < 20) . The areas are 
calculated as a sum of triangular components within the 
image/source polygon, where each triangle has a common 
vertex at y c or x c (ie. the null geodesic). 

We need to first check that the calculated /^true is not 
significantly in error using a particular value of i V rav . This 



was achieved by numerically integrating Equation ( A7), and 



comparing with fMrae for 4 < -/V ray < 256. The relative error 
in jUtrue was found to be independent of the choice of iV ray , 
and was well below the percentage cut-off level selected for 
fiaBM at the same radii. 

For source positions which are near the lens axis, it 
is difficult to keep track of which solutions x± belong to 
which of the images, particularly when merging of images 
is taking place. In these cases, fitmc is a misnomer, as it 
can be sig nific antly higher than the total magnification from 
Equation (A7), as may be seen in Figure ^. /irbm is subject 
to a similar error, as the areas here are also being calculated 
on the basis of iVr ay rays. We can proceed under the assump- 
tion that /itruc is sufficiently accurate for comparison with 
Mrbm when the two images are separable. 

By randomly selecting 5000 source (RBM image) lo- 
cations, we calculate fMrae and ^rbm- Table [l] shows the 
dimensionless radius, Scut, within which A/iRBM/Mtruc > 
p per cent for a range of source radii (10 -5 < Ri < 0.01). 
For a given value of the relative error, these results are es- 
sentially independent of both the source radii and iV ra y in 
the range 4 < iV ray < 256. 

The choice of iV ray = 8 appears to be the best compro- 
mise between accuracy and speed for our later investigations 
of the RBM. The smaller N Yay , the fewer total deflection cal- 
culations which need to be made, yet we retain high accuracy 
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0.2 0.5 1 2 

Source position: y 

Figure 3. Magnification of an extended source (R B = 0.4), near 
a Schwarzschild lens. /4rbM (thick solid line) is the Ray Bun- 
dle method magnification based on a single image; fitrue (thin 
solid line) is the total magnification for b oth i mages and ^integral 
(dashed line) is the solution of Equation (A7). /^rbm is accurate 
to 10 per cent at y cu t = 1.17, as indicated by the vertical dotted 
line. 



Table 1. Position of source, y cu t, at which the relative error in 
A'RBM (compared to the true magnification) is first > p per cent, 
caused by neglecting flux from second image. x cu t is the cor- 
responding (bright) image position. Values in the table are for 
source radii 10 -5 < Ri < 0.01 and 4 < 7V ra y < 256. 



fcut 
^cut 





^^RBm/ Mtrue 


= p% 




1% 


2% 5% 


10% 


20% 


2.86 


2.26 1.61 


1.15 


0.70 


3.17 


2.64 2.09 


1.73 


1.41 



for hhbm in the weak lensing limit. With 8 rays, we have a 
symmetric image with a quadrupole component, so it is pos- 
sible to determine the distribution of image ellipticities. 
In Figure 0, we plot /irbm (thick solid line), / itruc 



(thin solid line) and /^integral, the solution of Equation ( A7) 
(dashed line), as functions of the source impact parameter 
(y — \y\)- The source/image bundle radius in this case was 
R s — Ri — 0.4. We have used a larger radius than would be 
practical for the RBM in order to show in more detail what 
the high magnification behaviour is like near y = 0. As the 
source is comparable in size to the Einstein radius, the cal- 
culated magnifications near y ~ 0.4 are noisy. In this case, 
a circular source should produce a pair of highly distorted 
arcs, but this is not well represented by the use of trian- 
gular components within the ray bundle (ie. parity changes 
can occur in the individual triangles which make up the 
bundle). The vertical dotted line shows where the relative 
error in /irbm is accurate to 10 per cent, at y cut = 1.17 
or x cu t — 1.74 (note that these values are higher than in 
Table hj due to the larger source radius). 

This test has shown that the Ray Bundle method can 
give a highly accurate value of the magnification (relative 
error better than 5 per cent) in the weak lensing limit. How- 
ever, we have modified the Ray Bundle method slightly by 
requiring that the image position agrees with the brighter 
image of the corresponding circular source. In a true ap- 




0.5 1 1.5 

Impact parameter: x 



Figure 4. Magnification as a function of the impact parameter 
to a Schwarzschild lens. The thick solid line is the magnification 
for the RBM, the dashed line is for a point source with one image 
at the same location as the RBM image, and the thin solid line is 
the 'corrected' RBM magnification: fi' = Mrbm + 1- The Einstein 
radius corresponds to x = 1, the image radius is Ri = 10 — 4 and 
AW = 8. 



plication of the RBM, we randomly select image positions 
without the a priori knowledge of source locations. 

Figure ^ shows the RBM magnification for bundle po- 
sitions near, and within, the critical curve of a single 
Schwarzschild lens (thick solid line). For comparison, the 
total magnification of a point source with one image at the 
same location as the RBM image is shown as the dashed 
line. Bundle positions within x — 0.5 produce magnifica- 
tions /xrbm <C 1- For these bundles we actually selected 
the fainter of the images, as opposed to the case consid- 
ered previously, where we purposefully selected the bundle 
corresponding to the brighter image. Solutions of the lens 
equation, Equation (A3), which produce an image near the 
lens (x- < 1) correspond to total magnifications near (i = 1, 
as the source is far from the caustic point (y ~ —1/x- > 1 
This is demonstrated with the thin solid line in Figure 
where we plot ray bundles with \x < 1 within the Einstein 
radius as (/ = fi + 1. 

It is clear that a small strip of image locations near the 
critical curve is responsible for the high magnification region 
of the MPH, for which the Ray Bundle method is not well 
suited. 

Images well within the Einstein radius correspond to 
sources far from the lens axis. Such images are the faint 
images described in Section 



3.1 



(at Xfaint), and so there will 
be a second image at alright which contributes the majority 
of the flux. Although it is not possible to solve for alright 
given Xfaint, if the entire image plane is well sampled with 
ray bundles, then we can expect that another bundle will 
pass through Xbright and the source magnification will be 
calculated on the basis of this second bundle only. 

Two restrictions are now imposed on the Ray Bundle 
method for its later application to ensembles of lenses, which 
serve to complete the definition of the method. Firstly, im- 
age positions within the Einstein radius, or equivalently, any 
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magnifications which are calculated to be fi < 1 are dis- 
carded. Secondly, after selecting a relative error for ^rbm, 
we do not include images which fall within x C ut Einstein 
radii of any given lens. 

Having shown that in the 'weak' lensing limit (|x c | ;J 2), 
we can be sure of calculating the magnification to within 5 
per cent (or better) using TV ra y > 4 we can now proceed to 
a statistical comparison between the Ray Shooting and Ray 
Bundle methods. 



3.3 Comparison of magnification probability 
histograms 

Next we compare the MPH obtained with the RBM and the 
RSM. A number of subtle differences exist between the two 
methods, even when applied to the same lens model. The 
source size investigated is limited by the number of pixels 
in the source plane for the RSM. For a grid of TVpi x x TV p i x 
covering a square region 2t/ max x 2i/ max , the source 'radius' 

is Rs = ymax/iVpix. 

We choose the RSM sources to be squares with a side- 
length equal to the diameter of the circular RBM image bun- 
dles. Since the magnification decreases as the source area is 
increased, we expect each RSM source to have a system- 
atically smaller magnification than the corresponding RBM 
image. It is possible to use a much finer resolution grid for 
the ray shooting, and then integrate over a larger source 
size, but we have elected not to do this. This decision was 
based on a comparison of the computation time: for our im- 
plementation of the RSM, a grid of 1000 x 1000 pixels, with 
an average of TV = 250 rays per pixels took approximately 
eight hours of computation time. An equivalent number of 
images (where TV ray = 8) is completed in 1 minute with the 
RBM0 

Due to the distortion of rays near the boundary, as 
shown in Figure we actually shoot rays with the RSM 
through a larger angular region in the image plane. This 
prevents us from including source pixels which are not well 
sampled by rays. 

The accuracy of the RSM magnifications depends on 
the average number of rays collected in each pixel, TV. Since 
we are not implementing a (fast) hierarchical method, the 
time required to obtain the magnification distribution is pro- 
portional to the total number of rays, TVrsm = TVTVp ix . The 
RSM does not fully sample the highest magnification region 
(/i « //max) due to the regular placement of sources on a 
grid. 

With the RBM image size fixed by Ri — Rs, we are 
left to choose the number of rays which make up the ray 
bundle, TV ray , and the number of images. Ideally we want 
the images to completely cover the image plane (which again 
has a slightly larger angular size than the source plane due 
to the deflection of light rays near the boundaries) which 
requires 



Source Plane 



f 



irR? 



(2y n 



(14) 



t More computationally efficient versions of the RSM are avail- 
able, which can produce the same level of resolution in a 
time comparable to that of the RBM (J. Wambsganss, private 
communication) . 




Figure 5. A rectangular boundary in the image plane is dis- 
torted in the source plane. For the RSM, we only select sources 
in a regular grid which lies completely inside the distorted bound- 
ary, as shown by the dashed square in the right-hand panel. The 
Einstein radius is shown for scale. 



Table 2. Parameters for Magnification probability histogram 
comparison between Ray Bundle and Ray Shooting methods. 
/i m i n and //max are the minimum and maximum magnifica- 
tions calculated for each method, N( h B \ is the number of bun- 
dles/sources used for the RBM/RSM histogram, R(b, s ) ls the ra- 
dius of a circular bundle for the RBM and the half pixel size 



(see Figure p 


) and 


a distribution of five Schwarzschild lenses (see 


Figure H). In both 


case, TV = 250 


rays for the RSM and TV ray = 8 


for the RBM 












Method 


1 ' lens 


(Pmin - 1) 


^max 




Radius 


RSM 


1 


-0.03 


86.0 


10 6 


0.01 


RBM 


1 


2.4xl0~ s 


25.6 


9.8x10 s 


0.01 


RSM 


5 


-0.094 


37.7 


10 6 


0.01 


RBM 


5 


0.0 


218.1 


9.3 x10 s 


0.01 



where / is (approximately) the fraction of the source plane 
covered by sources. For the RBM then, the total number of 



rays required is TVrbm = TVi n 



x TV ri 



The comparative computational speed and result- 
ing magnification accuracy may be obtained by requiring 
TVrbm = TVrsm- Setting / = 1 and Ri = R B , we have 



TV = -TV ri 



We have already seen that TV ri 



(15) 

8 is a suitable choice, 
which gives TV ~ 10.2. This corresponds to a relative er- 
ror TV -3 / 4 ~ 17 per cent. Although we can relax the con- 
straint on the fraction of the source plane covered somewhat, 
and still have a well sampled MPH with the Ray Bundle 
method, we do not have this flexibility with the grid based 
Ray Shooting method. In addition, as we decrease the source 
size, the number of pixels required for the RSM increases, 
and a higher density of rays is necessary. 

Figure ^| shows the Magnification Probability his- 
tograms obtained for a single Schwarzschild lens using the 
Ray Shooting (thin line) and Ray Bundle (thick line) meth- 
ods. The parameters for each method are listed in Table ^. 
The vertical axis of this (and later histograms) is the nor- 
malised number of bundles/sources in each magnification 
bin, N(fi), which is equivalent to the definition of p(/ii)A/i 
in Equation (Q). 

A cut-off was imposed on image locations for the RBM 
at s cut = 1.01 Einstein radii. The two distributions are 
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Figure 6. Comparative magnification probability histograms 
for the Ray Bundle method (thick solid line) and the Ray Shoot- 
ing method (thin solid line) for a single Schwarzschild lens. The 
dashed line shows the expected high magnification power law 
slope N(fj,) oc pL~ 2 . See Table for the model parameters. 

qualitatively very similar, when the various caveats de- 
scribed above are considered. Imposing a larger value of 
Xcut serves to reduce the maximum magnification with the 
RBM. The poor sampling in the highest \i bins for the 
RSM is clearly demonstrated. The dashed line shows the 
expected fi~ power law slope of the MPH at large \i (see 
Appendix |a|) , and both distributions have this approximate 
form (although for the RSM the statistical significance of 
the histogram bins with /i > 10 is low). 

As expected, the RSM produces magnifications which 
are fi < 1 (when /x cmpt y > 1 is expected) due to numerical 
effects. The sample mean and variance of the two distribu- 
tions are (fi) — 0.98 and a 2 . = 0.03 for Ray Shooting, and 
(n) = 1.02 and = 0.08 for the Ray Bundle method^} 
The Ray Shooting method provides higher accuracy at high 
magnifications, but at magnifications fj, ~ 1, the Ray Bundle 
method is more accurate even though flux from additional 
images is neglected. 

One aspect of the MPH we have not yet discussed has 
to do with the weighting we apply to each ray bundle. For 
the RSM, every source is a pixel with the same area. For the 
RBM, the initial ray bundles (images) have the same area, 
but the resulting sources must have different areas (by the 
definition of a magnification). It is sufficient to weight each 
ray bundle in the MPH by the area of the resulting source. 
Figure Q shows the RBM (thin line) with no weighting com- 
pared with the correct area weighting (thick line). 
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Figure 7. Dependence of the Ray Bundle method magnifica- 
tion probability histogram on the weighting. With no weighting 
(thin solid line) the MPH is incorrect at high magnification, where 
the source sizes are much smaller than in the weak lensing limit, 
while weighting by source area (thick solid line) gives the expected 
power law slope N(fi) oc u~ 2 at fj, > 10 (dashed line). See Table ti 
for the RBM model parameters. 

3.4 More complex lens distributions 

For an ensemble of N lenses in the lens plane, each with 
mass Mj, the total deflection angle generalises to 

where the (^ — £j) are the impact parameters to each lens. 

Consider the case where lenses are restricted to lie 
within a rectangular region^]. Light rays passing through one 
of the corners of the shooting region will necessarily be de- 
flected inwards by the mass distribution. This is appropri- 
ate for an isolated configuration of lenses, such as in studies 
of the microlensing effect of many stars which make up a 
galaxy (where the contribution to the deflection by an ex- 
ternal mass distribution may be modelled by adding a shear 
term to the lens equation). However, for an investigation 
of the lensing due to large scale structure, where the mass 
distribution is assumed to be continuous and homogeneous 
in all directions about a ray, we may introduce an artificial 
shear on rays near the shooting boundary. 

For the RBM, we choose to calculate each of the deflec- 
tion angles explicitly, with an increase in the computational 
time over an equivalent hierarchical method. By making the 
direct calculation of the deflection, we are free to choose 
the geometry of the region within which we include lenses. 
The most natural choice for a distribution which is homo- 



+ We reiterate that in the definition of the RB M. we do not in- 
clude images with fj, < 1, see end of Section 3.2). 



3 This is the case in common implementations of Ray Shooting as 
i t allows for the easy implementation of fast, hierarchical methods 
(Wambsganss 1990[). 
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Figure 8. Comparative magnification histograms for the Ray 
Bundle method (thick solid line) and the Ray Shooting method 
(thin solid line) for a distribution of 5 equal-mass Schwarzschild 
lenses. The dashed line shows the expected high magnification 
power law slope N(fi) <x /J -2 . See Table ti for the model parame- 
ters. 



geneous and isotropic beyond some length scale Rh is to 
include lenses within a circular region around each ray out 
to the radius 7?i ens = Rh- For the isolated lens geometries 
we examine here, we need only set Ri ens to encompass the 
RSM image plane. A discussion of the appropriate choice for 
Rh in cosmological lensing scenarios is reserved for Fluke et 
al. (in preparation). 

As a final test here, we now consider a random distri- 
bution of iViens = 5 lenses (with the same lens positions for 



both th e RSM and RDM). The resulting MPII is shown in 



Figure g|and variouo paramotora in Table For both moth 

ods we use a total of 10 6 bundles/sources, but discarding 
bundles within a; cut = 1.01 Einstein radii for the RBM re- 
duces this to 9.3 x 10 bundles for the histogram. The sample 
mean and variance are (fj) = 1.09, a 2 , = 2.09 for the RBM 
and (fj) = 0.98, a% = 0.23 for the RSM. 

Using an ensemble of lenses introduces a new length 
scale to the problem, so that sub-structure at intermediate 
magnifications (3 ~ fi ~ 30) is seen in the MPH using the 
RSM (but not with the RBM). The 'bump' in the MPH oc- 
curs for a planar d istribution of lense s, and has been st udied 



both n umerically ( Flauch et al. 1992 ) and analytically ( Kof- 
man et al. 1997), and is believed to be due to the caustic 



patterns of pairs of point lenses. If this is a caustic-induced 
feature, ie. a region of high magnification, then the RBM 
will not provide the 'correct' magnification. We feel the fea- 
ture may be in part due to the low resolution with which 
the complex caustic structure was mapped with the RSM 
on a regular grid (1000 x 1000 pixels), however a discussion 
of this effect is beyond the scope of this paper. 



4 CONCLUSIONS 

The Ray Bundle method provides a computationally fast, 
accurate and flexible alternative to the Ray Shooting 
method for studies of the weak gravitational lensing limit. A 
wide variety of lens models are easily incorporated - we have 
considered here only the case of the Schwarzschild lens, but 
changing to a different model involves modifying only «(£) 
in the GLE. One alternative to ray based methods requires 
solving for the scalar lensing potential, but it may not always 
be possible to find an analytic solution for all lens models. 
The RBM also allows us to avoid artificial shear introduced 
by grid based methods for light rays which pass near the 
shooting boundary, when a small portion of an otherwise 
homogeneous and isotropic distribution is used. 

The point source limit may be approached as values of 
Ri may be selected without the restriction of introducing a 
finer source grid, and a corresponding increase in the density 
of rays required (provided we can relax the constraint that 
the source plane must be completely covered by sources with 
the RBM). 

The RBM should only be applied with caution to strong 
lensing scenarios, where the RSM is far superior. As only one 
image is followed to the source, there will be an error in the 
total magnification when the image is near a critical curve. 

The RBM is, however, particularly well suited to prob- 
lems were we want to investigate in detail individual lines-of- 
sight for various lens geometries and models. An important 
advantage of the RBM is that we can associate a particu- 
lar image position and shape with the corresponding source 
position and shape. This provides us with the opportunity 
of following the development of the shape of a ray bun- 
dle through a sequence of lens planes, as used in models 
of cosmol ogica l lensing (for example, Wambsganss, Cen & 
Ostriker (p98h). 



An important application of weak lensing is to deter- 
mine the effect of small changes in the magnification of 
standard candle sources (such as Type la Su pernovae) on 
the derived values of cosmological parameters ( [Wambsgans;; 



et al. 1997). The high accuracy of the RBM in the weak 



lensing limit makes it a value tool for such studies. 
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If the histogram bins are equally spaced logarithmically, 
in 1^2 is constant, and so we have an additional constraint 
that the MPH must have a power law slope —2 for /i > 10. 

For an extended source, the total magnification is the 
integral of /i p over the source, weighted by the intensity 
profile, T(y). For a circular source with dimensionless ra- 
dius Rs = Rs/rjp, and a uniform intensity profile we have 



formulation) 



(eg. Pei (1993), or Schneider et al. (1992) for an alternative 



y+Rs 
dt- 



*Ri 4_k s | y/R2-{y-t)*y/&+W 2 



Rl 



(A7) 



Equation (A7) approaches the point source solution as 
y — > oo or R s — » 0, and the maximum magnification is 



/Xe,max — 



R s 



(A8) 



at y = 0. 



This paper has been produced using the Royal Astronomical 
Society/Blackwell Science KLgX style file. 



APPENDIX A: THE SCHWARZSCHILD LENS 

The simplest lens model is the po int-m ass or Schwarzschild 



lens ( for ex ample Schneider et al. (1992), Narayan & Bartel- 



man ( 199!: )) , for which the deflection angle due to a mass, 
M, is 

"«) = ^J2^ ( Al ) 

The dimensionless lens equation for a point source with this 
model is 

y = x-l/x (A2) 

so that there are two images (one located on either side of 
the lens, and co-linear with the lens) at 

x± = \{y±yJv* + ±) (A3) 
with corresponding magnifications 

«-i(-£i_ ± lY (A4, 

2 \yVy + 4 J 

The total magnification is the sum of the absolute values of 
the individual image magnifications: n P = ^t+| + |/i- |. When 
the source is far from the lens axis (y > 1), one of the images 
(say, X-) will be significantly demagnified (^t_ <C 1). The 
total magnification is then fi p ~ 

It is possible to derive an analytic form for the mag- 
nification probabili ty, p(p,,z), for the Sch warzschild lens for 



large values of fi p ( Schneider et al. 1992), which is reason- 
ably generic for most lens models ( [Peacock 1982] ): 

p(p) oc n~ 3 . (A5) 

On integrating to form the MPH we have 

1 {, A*?' 



p( M )A M = / p( M )d M = — 1 - £^ ) . (A6) 
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